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While phylogeographic data provide valuable information to inform conservation plans, there are com- 
paratively few Odonata phylogeographic studies. This lack of research is partially due to a lack of 
independent DNA markers with appropriate levels of polymorphism that PCR-amplify in a range of 
species. We followed an exon-primed, intron-crossing (EPIC) PCR strategy to develop five new, poly- 
morphic nuclear DNA sequence loci (six distinct DNA fragments) for the southern damselfly Coenagrion 
mercuriale. These markers were: cell division cycle 5 protein (CDC5), arginine methyltransferase 
(PRMT), acetylglucosaminyl-transferase (AgT), myosin light chain (MLC) and phosphoglucose iso- 
merase (PGI). Between three and five of these new markers could be PCR-amplified in five other species 
from the genus Coenagrion; one locus (PRMT) can be used in 26 other species of odonates that we exam- 
ined, including three species of Anisoptera belonging to the genus Onychogomphus. These new nuclear 
genetic markers will be useful for phylogeographic studies in a range of odonate species, but also for phy- 
logenetic studies, providing a particularly useful complement to the existing mitochondrial and nuclear 
loci. 
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Introduction 


Phylogeographic studies yield insights into the suite of historic processes that have shaped the 
contemporary distribution of genetic diversity (Avise, 2000). Such information assists effec- 
tive conservation management by providing fundamental data about evolutionarily significant 
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units (ESUs), thereby allowing prioritization of areas/populations that represent high value for 
conservation (Moritz, 2002). 

Regions of mitochondrial DNA (mtDNA) are still the most widely used molecular-genetic 
markers to assess phylogeographic structure in animal populations (Avise, 2000, 2009; Ballard 
& Whitlock, 2004). The popularity of mtDNA for phylogeographic studies is largely due to (1) 
its higher rate of mutation than in most nuclear genes; (2) a predominantly maternal inheritance 
that generally excludes recombination as a source of diversity and ensures homoplasmy; and 
(3) a conserved and tightly packed gene content that allows development of universal primers 
for PCR (e.g. Simon et al., 1994). While the mutation rate varies among mitochondrial genes, 
resulting in distinct phylogenetic resolution power, the use of mitochondrial genetic markers can 
be problematic to interpret. For example, the phylogeographic history of an organelle (such as 
the mitochondrion) may differ from an organism’s history due to selective sweeps and intro- 
gression, or due to demographic attributes, such as population size, dispersal ability and/or 
sex-biased dispersal (Ballard & Whitlock, 2004; Galtier, Nabholz, Glémin, & Hurst, 2009; Hare, 
2001). Amplification of numts (nuclear mitochondrial DNA; Lopez, Yuhki, Masuda, Modi, & 
O’Brien, 1994), where whole or fragments of mitochondrial genes are incorporated into the 
nuclear genome, might constitute an additional problem by providing a heteroplasmic signal 
(in the absence of true heteroplasmy). Nevertheless, with careful analysis numts can be used to 
uncover complex evolutionary histories in the context of phylogeographic analysis (e.g. Miraldo, 
Hewitt, Dear, Paulo, & Emerson, 2012). Further problems with the use of mtDNA as a phy- 
logeographic marker may arise in studies of arthropods due to indirect selection on mtDNA 
arising from linkage disequilibrium with maternally inherited symbionts (Hurst & Jiggins, 
2005). Robust phylogenetic inference thus requires data from several independent genetic 
markers. In animals, such independent genetic loci must be derived from the nuclear genome 
as mitochondrial genes are physically linked and thus represent a single locus, except for 
some rare exceptions (e.g. Jiang, Barker, & Shao, 2013). Nonetheless, an apparently low 
rate of mutation found in most nuclear genes remains a major drawback for their use in 
phylogeography. 

Introns (non-coding regions of genes located between exons) exhibit higher levels of poly- 
morphism than exons (segments of genes that code for a peptide) and thus provide a potentially 
valuable source of nuclear genetic variation for phylogeographic studies. However, it is challeng- 
ing to design reliable PCR primers in polymorphic regions of the genome, particularly if the goal 
is to amplify the locus in several species. The use of exon-primed intron-crossing (EPIC) mark- 
ers resolves this issue by targeting conserved regions of the genome (i.e. the exons) to anchor 
primers that are used to PCR-amplify and sequence the adjacent introns (Lessa, 1992; Slade, 
Moritz, Heideman, & Hale, 1993). Indeed, EPIC markers have been successfully used in phylo- 
geographic studies of a variety of taxa (e.g. Lohse, Sharanowski, & Stone, 2010; Aurelle et al., 
2011; Jennings, Etter, & Ficarra, 2013). 

Phylogeographic analyses of odonates remain relatively scarce compared with other animals 
(but see Artiss, 2004; Bernard, Heiser, Hochkirch, & Schmitt, 2011; Biisse et al., 2012; Jordan, 
Simon, Foote, & Englund, 2005; Kiyoshi, 2008; Kiyoshi & Sota, 2006). Notably, these odonate 
studies continue to rely upon a limited set of genes: cytochrome c oxidase I (COD, cytochrome 
c oxidase II (COII), NADH dehydrogenase 1 (ND1) and 16S ribosomal DNA (16S), and the 
nuclear ribosomal RNA gene region encompassing 18S, ITS1, 5.8S and ITS2. The present 
reliance on using mtDNA markers seems likely to continue with the recent publication of some 
entire odonate mitochondrial genomes (e.g. Chen et al., 2014; Lee et al., 2009; Lin et al., 2010; 
Lorenzo-Carballa, Thompson, Cordero-Rivera, & Watts, 2013; Wang et al., 2014). While these 
data offer the prospect of identifying universal PCR primers for the amplification of additional 
mtDNA genes in odonates, they do not resolve the key issue that all genes located on the mtDNA 
represent a single locus. 
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Just seven nuclear DNA sequence markers (18S, ITS1, 5.8S, ITS2, 28S, H3 and EF-1a@) have 
been widely employed to quantify inter- and intra-specific relationships in odonates (Ballare & 
Ware, 2011). However, five of these markers (18S, ITS1, 5.8S, ITS2 and 28S) are physically 
linked to represent a single locus; moreover, the ribosomal coding sequences (5.8S, 18S and 
28S) and histone 3 (H3) are conserved and thus largely uninformative for intraspecific phylo- 
geographic studies (e.g. Cruickshank, 2002; Simon, Schierwater, & Hadrys, 2010). Elongation 
factor-1a@ (EF-1œ) has been used to resolve odonate phylogenies (e.g. Carle, Kjer, & May, 2008; 
Jordan, Simon, & Polhemus, 2003) but its suitability to recover phylogeographic patterns is 
unclear (Simon et al., 2010). Essentially this leaves the ribosomal internal transcribed spacers 
(ITS1 and ITS2) that present intraspecific polymorphism in odonates (e.g. Dumont, Vanfleteren, 
De Jonckheere, & Weekers, 2005; Weekers, De Jonckheere, & Dumont 2001), thus there is a need 
to develop additional loci (1) to achieve statistical confidence and (2) because the ITS is a multi- 
copy gene that can exhibit substantial intragenomic variability which obscures phylogeographic 
patterns (Harris & Crandall, 2000). 

At present, only Lee and Lin (Molecular Ecology Resources Primer Development Consortium 
et al., 2012) have attempted to develop new nuclear DNA sequence markers in odonates. Eight 
EPIC loci were developed for Euphaea spp. (Euphaeidae) but their utility in other odonates was 
not tested. Thus, Odonata phylogeographic studies remain hampered by the lack of a set of inde- 
pendent, polymorphic nuclear DNA markers that are useful across a range of taxa. To redress this 
deficiency and to improve our understanding of odonate phylogeographic patterns using a mul- 
tilocus approach, we developed five new nuclear DNA sequence loci in the southern damselfly 
Coenagrion mercuriale. These loci can be used in other zygopterans and in some anisopter- 
ans (Onychogomphus spp.), representing a wide taxonomic diversity encompassing regions in 
Europe, North Africa and Central America. 


Material and methods 


Target species and cross-amplification 


Coenagrion mercuriale is one of the 19 Western Palaeartic species of the genus Coenagrion, 
and belongs to Coenagrionidae, the largest family of zygopterans with almost 1100 species 
distributed on all continents except Antarctica (Dijkstra & Kalkman, 2012). Coenagrion mer- 
curiale’s distribution extends from the UK and Germany through France, Iberia, Italy and into 
North Africa (Dijkstra & Lewington, 2006). This species is classified in the IUCN red lists as 
Near Threatened in Europe (Kalkman et al., 2010) and as Endangered in North Africa (Samraoui 
et al., 2010). We sequenced 10 individuals of C. mercuriale (five from Europe and five from 
North Africa) to evaluate polymorphism in the newly developed markers. 

Exploratory cross-amplification tests using the PCR primers for the new loci developed for 
C. mercuriale were performed in 26 other species of zygopterans and anisopterans (Table 1). 
Within the coenagrionids, we selected five species of Coenagrion (C. caerulescens, C. orna- 
tum, C. puella, C. pulchellum and C. scitulum), one species of Pseudagrion (P. sublacteum) 
and six species of Ischnura (I. hastata, I. ramburi, I. cervula, I. denticollis, I. kellicotii and 
I. posita). Representatives of three other zygopteran families were Protoneura capillaris, P. 
viridis, P. sanguinipes, Neoneura amelia, N. maria and Microneura caligata (Protoneuridae; 
recently suggested to be part of Coenagrionidae s. lat., Carle et al., 2008; Dijkstra et al., 2014), 
Hypolestes clara and H. trinitatis (Hypolestidae) and Calopteryx haemorrhoidalis, C. xanthos- 
toma, C. virgo and C. splendens (Calopterygidae). Cross-amplification tests were also performed 
in three species of anisoptera from the genus Onychogomphus (O. costae, O. forcipatus and 
O. uncatus). 


Table 1. 


Genbank accession numbers. 


Amplification success of six new EPIC nuclear loci in 27 odonate species displayed against a schematic representation of their putative evolutionary 
relationships (derived from Dijkstra et al., 2014). Data indicate the number of specimens sequenced, the primer pairs to PCR-amplify gene fragments and the 






































































































































cocs MLC PRMT Agr PGES PGE1219 
Genus N Primers Genbank N Primers Genbank N Primers Genbank N Primers Genbank N__ Primers Genbank N Primers ank 
mercuriale 10 O FèR2  KM248413-KM248432? | 10 Flu-Ricm_  KM248443-KM248462 | 10 _F4-R3 10 FRA  KM248393-KM248412 | 10 8F-382R_  KM248463-KM248482 | 10 1I219F-1455R__KM248483-KM248502 
ornate 1 Flu-Ricm KM276626 1 FAR KM276639 1 F2-R4 KM276647 1 8F.382R KM276651 1 1219F-1455R KM276638, 
PORRA paik L F2-R2 KM276596 L__FluRlem KM276623 1 FAR KM276640 1 PR KM276644 1 BF-382R KM276652 1 1219F-1455R. KM276637 
ananin 1 PR KM276598 L__FluRlem KM276625 1 FAR KM276641 Lore KM276646 1 8F-382R KM276650 1 1219F-1455R KM276635 
seitulum 1 Flu-Rlem KM276627 1 FAR3 KM276643 1 8F-382R KM276653 
aende L F2-R2 KM276597 1 Flu-Ricm KM276624 1 FR KM276642 1 PR KM276645 1 BF-382R KM276654 1 1219F-1455R. KM276636 
Pseudagrion sublacteum 1 O FR KM276649 
hastata 9 FARI KM276587-KM276595 16 ARS KM276605-KM276620 
cervula 1 FaR3 KM276622 
Ischnura denticollis 1O FR KM276603 
kellicoit 1 FR KM276621 
posita 1 O FR KM276604 
capillaris 1 F2-R3 KM276631 
viridis 1 PR KM276632 
sanguinipes 1 PR KM276634 
Microneura caligata 1 F2-R3 KM276633 
T amelia 2 F}R3  KM276629-KM276630 
maria 1 FR KM276628 
a pan 1 BB KM276601 1 BF-382R KM276599 
trinitatis 1 Farid KM276648 1 mR KM276602 1 KM276600 
haemorrhoidalis 3 F}R3  KM276570-KM276572 
PATNA splendens 5 F2-R3 KM276582-KM276586 
virgo 6 F}R3  KM276573-KM7 
xanthostoma 3 F}R3  KM276579-KM 
uncas 1 FR KM222659 
Onychogomphus | forcipatus 1 OPR KM222656 
costae 1 PR KM222662 
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Marker selection and development 


As our goal was to identify putative single copy nuclear DNA sequence markers that were suf- 
ficiently polymorphic for phylogeographic studies we adopted an exon-primed, intron-crossing 
(EPIC) PCR strategy (Lessa, 1992). From an arthropod phylogenomic analysis of 62 single-copy 
nuclear protein-coding genes (Regier et al., 2010), which included the odonates [schnura verti- 
calis and Libellula lydia, we first selected the genes considered as “sufficiently variable (...) to 
be a candidate for testing across Lepidoptera” (Regier, 2007, pp. 15-33) and then identified those 
genes with introns in other insects using the BLAST (Altschul, Gish, Miller, Myers, & Lipman, 
1990) tool in flybase (McQuilton, St Pierre, Thurmond, & FlyBase Consortium, 2012). We then 
downloaded sequence data from I. verticalis and L. lydia for this reduced set of nuclear genes 
from Genbank (http://www.ncbi.nlm.nih.gov/genbank). Due to apparently high polymorphism 
(Lee & Lin, 2012) we also downloaded the Euphaea sp. sequence data for myosin light chain 
(MLC), arrestin (ARR) and ferritin (FER) sequences from Genbank. 

To obtain sequence data for primer design, we used a BLAST search (Altschul et al., 1990), 
performed in Geneious v.6.1.5 (http://www.geneious.com), to identify which of the candidate 
genes were present in a C. mercuriale transcriptome (derived from 454 sequencing of RNA 
expressed in larval thoracic muscle tissue; P.C. Watts and D.J. Thompson, unpublished). The 
putative locations of introns in the C. mercuriale transcriptome sequences were identified by 
comparison against homologous sequences from other insects. Primers around introns were 
designed manually in Bioedit v.5.0.9 (Hall, 1999). 

Additionally, phosphoglucose isomerase gene (PGI) was examined as a candidate gene for 
PCR-optimization, and polymorphism evaluation in our target species using Coenagrion sp. spe- 
cific primers (Fox, 2010), as it is a candidate locus for dispersal ability in insects (Wheat, Watt, 
Pollock, & Schulte, 2006). For comparison, we sequenced two widely used nuclear sequence 
markers: a fragment of 28S gene (Hasegawa & Kasuya, 2006) and a fragment of H3 (Ogden & 
Whiting, 2003). 


Molecular biology and data analysis 


Total genomic DNA was isolated from leg samples using a high-salt protocol (Sambrook & 
Russell, 2001). PCRs were carried out in 10 pl reaction volumes containing 5 u1 of ReddyMix™ 
PCR Master Mix (Thermo Scientific, Waltham, Massachusetts, USA) [75 mM Tris-HCl (pH 8.8), 
20mM (NH4)2SOq, 1.5mM MgCh, 0.01% (v/v) Tween® 20, 0.2mM each dNTP and 0.625 
units ThermoPrime Taq DNA Polymerase], 1.5 mM MgCl, 0.2 pmol of each primer, and ~ 10 ng 
of genomic DNA. Thermal cycling conditions consisted of 5 min at 94°C, followed by 40 cycles 
of 30s at 92°C, 30 s at TC and 45 s at 72°C, and then 10 min at 72°C, where 7, is the marker- 
specific annealing temperature (Table 2). 

Following PCR, the remaining primers were removed by incubation with 0.75 units Antarctic 
Phosphatase (New England Biolabs, Ipswich, Massachusetts, USA) and 0.6 units Exonuclease I 
(E. coli) (New England Biolabs) at 37°C for 75 min, followed by 15 min at 80°C [in 7 ul final 
reaction volumes containing 1X Exonuclease I Reaction Buffer (New England Biolabs)]. PCR 
products were sequenced in both directions using BigDye v.3.1 chemistry (Applied Biosystems, 
Foster City, California, USA), following the standard protocol but with 10 ul reaction volumes. 
Purified sequence products were run on an ABI3100x1 (Applied Biosystems). 

Chromatograms were checked for quality using Geneious v.6.1.5 (http://www.geneious.com) 
and nucleotide ambiguities with similar peak heights were considered to be heterozygous posi- 
tions and were recoded with IUPAC ambiguity codes. Sequence alignment was conducted using 
MUSCLE (Edgar, 2004) (implemented in Geneious) with 10 as the maximum number of iter- 
ations, and alignments were then corrected manually if necessary. Haplotype reconstructions 
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Table 2. Details of PCR primers and thermal cycling conditions for successful amplification of regions of nuclear 
genes in odonates. T, represents the marker-specific annealing temperature. 





Gene (abbreviation) Primer Primer sequence 5’—>3’ Ta CC) 
Arginine methyltransferase (PRMT) ARG_F2 TGC CGC CAA GGC TGG AGC ATC 50-53 
ARG_F3 CCG GAA CTC TAT GTA CCA CAA C 
ARG_F4 TCG ACT CGT ATG CGC ATT TCG G 
ARG_R3 TGC CAC CTT CCT AAT AGA GCT C 
Phosphoglucose isomerase (PGI) Cp pgi 8 F GTC GTG CAA AAC TCA CAG AAG AGG 56 
Cp pgi 382 R TTC CTT CAT GTG ACT CAG AAC TGC 
Cp pgi 1219 F CTG CTG ACT TCA TAG CCC CTG TAA 56 
Cp pgi 1455R GGC CCC WAG AGT AAA AGG TGT GAC 
Myosin light chain (MLC) Myo_Flu ACT TCA CCC AAC TGC TCAC 47 
Myo_R1cm CAT CAT CGA ATG ACT TGA 
Myo_F2d CAC CCA ACT GCT CAC MCT TT 50 
Myo_Rld CAT CRT CGA ATG ACT TGA 


Cell division cycle 5 protein (CDC5) CDC5_F4 GCT CGT TGG TAY GAA TGG TTA GAT CC 45 
CDC5_R1 GCT GGC TTT GTT TCT GGA TTA GGA TC 
CDC5_F2 GCA AAA TTG ATG CCT ACM CAA TGG AG 45 
CDC5_R2 GGA TTA GGA TCA ATT TCT CCA GGC TTC 
Acetylglucosaminy]-transferase (AgT) AgT_F2 TGC GTG CTC TTA ACC TTA GC 53 
AgT_R4 TAG GCT GAA GCT CAA TTG CTC 


were performed using PHASE v.2.1 (Stephens & Donnelly, 2003; Stephens, Smith, & Donnelly, 
2001) implemented in DnaSP v.5 (Librado & Rozas, 2009). Unrooted haplotype networks were 
constructed using statistical parsimony (Templeton, Crandall, & Sing, 1992), to provide a visual 
comparison of the level of resolution of each genetic marker, using TCS v.1.21 (Clement, Posada, 
& Crandall, 2000), with indels treated as a fifth character state. 

We calculated standard summary statistics of diversity (S, number of segregating sites; 7, 
nucleotide diversity; h, number of haplotypes; Hd, haplotype diversity), the minimum number 
of recombination events (Ry) (Hudson & Kaplan, 1985) and respective statistical confidence 
based on a coalescent simulation algorithm, and tested for deviation from neutral predictions 
using Tajima’s test (D) (Tajima, 1989) using the phased data in DNAsp v.5.10.01 (Librado & 
Rozas, 2009). The phi-test (Bruen, Philippe, & Bryant, 2006) was performed to assess statistical 
evidence for recombination in SplitsTree v.4.13.1 (Huson & Bryant, 2006). 


Results 


Sixteen of the putative variable genes from Regier (2007) had one or more confirmed introns 
in insects. We found six of these genes [EF-1qa, cell division cycle 5 protein (CDC5), arginine 
methyltransferase (PRMT), acetylglucosaminyl-transferase (AgT), alanyl-tRNA synthetase, and 
transmembrane protein] in the C. mercuriale transcriptome, as well as sequences of the myosin 
light chain (MLC) and phosphoglucose isomerase (PGI). Sixty-five primers were designed for 
regions within these eight genes. We achieved successful PCR-amplification from C. mercuriale 
genomic DNA in five genes: CDC5, PRMT, AgT, MLC and PGI (Tables 1, 2,3). We sequenced 
the longest fragments that could be amplified for each gene, although for the PGI locus two 
non-overlapping fragments (PGI-8, PGI-1219) were sequenced (see Tables 1, 2, Figure 1). 

The five new nuclear loci sequenced present distinct structural features regarding fragment 
length, presence of introns and presence of indels (Table 3, Figure 1). The new markers are capa- 
ble of amplifying some 3200 bp in C. mercuriale, representing fragments ranging from ~200 bp 
(CDCS) to ~750 bp (PRMT and PGI-1219) (Table 3, Figure 1). In C. mercuriale, the fragment 
length varies in four loci due to the presence of indels, with the longest indel being 6 bp long in 


Table 3. 


Summary diversity statistics for eight nuclear sequence markers in 10 specimens of Coenagrion mercuriale representing individuals from Europe (n = 5) and North Africa 





(n= 5). 

FW primer RV Length Indel no./ Intron no./ 
Gene pair primer (bp) length (bp) length (bp) S T h Hd Rm w D 
2884 28sf 28sr 830 0 0 0 0 1-1 0 NA NA NA 
H38 HexAF HexAR 284 0 0 2 0.0038 2-2 0.556 0 1.000 ns 1.8443 ns 
CDC5 CDC5_F2 CDCS5_R2 208 0 1 (88) 5 0.0049 6 0.795 0 1.000 ns —0.3060 ns 
MLC Myo_Flu Myo_Ricm 319 1 (1) 1 (257) 16 0.0170 8-8 0.821 0 1.000 ns 0.8076 ns 
PRMT ARG_F4 ARG_R3 759 1(1) 3 (97, 78, 89) 17 0.0090 8-9 0.874 0 1.000 ns 1.3426 ns 
AgT AgT_F2 AgT_R4 603 3, 1,4) 1 (518) 22 0.0134 8-10 0.884 2 0.229 ns 1.2025 ns 
PGI Cp pgi 1219F Cp pgi 1455R 746 3 (6, 1, 2) 1 (>612) 17 0.0091 7-7 0.832 0 1.000 ns 1.4732 ns 
PGI Cp pgi 8F Cp pgi 382R 619 0 2 (145, 171) 8 0.0055 6-6 0.758 0 1.000 ns 1.3260 ns 


Length, length of sequence in the final alignment (bp); Indel no., the number of insertion/deletions polymorphisms and respective length; Intron no. the number of introns and respective length(s); S, number of 
segregating sites; 2, nucleotide diversity; h, number of haplotypes (calculated by DNAsp and TCS respectively); Hd, haplotype diversity; Ry, minimum number of recombination events; w, phi-test; D, Tajima’s 
test. ns, not significant; *p < 0.05. 
AHasegawa & Kasuya (2006) 
BOgden & Whiting (2003) 
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Figure 1. Diagrammatic representation of the six fragments of the five new, nuclear loci to indicate the locations of 
primer annealing sites relative to the exons and introns. The coding sequence (solid black line) represents C. mercuriale 
transcriptome sequence, while exons (dark grey) and introns (light grey) represent new sequence data. The first 10 lines 
represent individuals of C. mercuriale, whereas each of the following lines represents successful amplification of other 
species sequences during cross-amplification tests (see Table 1 for details on species). 


PGI-1219 (Table 3). As we followed an EPIC strategy, all new markers encompass between one 
and three introns, which represent the location of most polymorphisms, with the smallest intron 
being 78 bp in length (in PRMT) and the longest intron more than 612 bp (in PGI-1219) (Table 3, 
Figure 1). 

All new markers exhibited instances of heterozygous individuals (in the C. mercuriale dataset) 
and presented a significant increase in polymorphism when compared with the commonly used 
nuclear fragments (28S and H3) (Table 3, Figure 2); thus the new loci provide good phylogeo- 
graphic resolution. Among the 10 European and North African C. mercuriale, the new markers 
yielded between two and five haplotypes within regions, which contrasts with the 28S and H3 
fragments where there were just one or two haplotypes in total (Figure 2). Moreover, the Euro- 
pean and North African samples are genetically distinct, and reasonably divergent, at all new 
markers, except CDC5 in which European and North African specimens share the most common 
(possibly ancestral) haplotype (Figure 2). 

No deviation from neutral predictions was apparent at any locus (i.e. Tajima’s D), although 
the sample size is small (Table 3). Recombination was detected at AgT, with a significant (p < 
0.05 after coalescence simulations) minimum number of recombination events of 2, although the 
phi-test did not uncover evidence of recombination in any locus (Table 3). 

Successful cross-amplification and sequence data were obtained for all new markers for four 
species of Coenagrion, using the same primers that proved optimal for C. mercuriale (Table 1) 
with C. scitulum failing to amplify two fragments (AgT and PGI-1219). Arginine methyltrans- 
ferase (PRMT) was successfully amplified in all 27 odonates species examined. The PCR primer 
pairs that amplify this locus in coenagrionids and Hypolestes species generates fragments of 
~760 bp and ~700bp in length, respectively, and which encompasses three introns (Table 3, 
Figure 1). In all other species, including the anisopterans, successful PCR amplification of PRMT 
could only be achieved with the primer pair that targets a smaller (~500 bp) fragment that encom- 
passes two introns (Table 3). Other loci also amplified fragments in non-coenagrionids, most 
notably MLC and PGI-8 in both Hypolestes species, reinforcing the potential for the use of these 
loci in other odonates (Table 1). 
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Figure 2. Haplotype networks (generated using TCS v.1.21) indicating diversity and relationships among Coenagrion 
mercuriale from Europe (white, n = 5) and North Africa (grey, n = 5) to illustrate the level of phylogeographic resolution 
yielded by eight nuclear DNA markers (see Tables 1, 2 and 3 for details of loci). Circle size is proportional to haplotype 
frequencies and black points represent missing or unsampled haplotypes. 








Discussion 


Phylogeographic studies of odonates remain relatively scarce albeit with increasing prominence 
in recent years (e.g. Artiss, 2004; Bernard et al., 2011; Büsse et al., 2012; Hayashi, Dobata, 
& Futahashi, 2005; Jordan, Simon, Foote, & Englund, 2005; Kiyoshi, 2008; Kiyoshi & Sota, 
2006). As phylogeographic analyses of odonates typically rely on mtDNA markers, they leave 
unexplored the type of information and associated advantages that can be derived from nuclear 
markers. But the crucial need for independent markers is highlighted by reports that phylogeo- 
graphic patterns based on mtDNA markers can be incongruent with those derived from nuclear 
markers (e.g. Hayashi et al., 2005). The absence of a set of independent DNA markers with 
appropriate levels of polymorphism has been an apparent constraint for making robust studies of 
intraspecific diversity in odonates. Using an exon-primed, intron-crossing (EPIC) PCR strategy 
we developed five new, polymorphic nuclear DNA sequence loci to improve the resolution of 
phylogeographic studies of coenagrionids and other odonates. 

Besides the key issue of marker independence, the use of nuclear DNA sequence markers 
allows the detection of hybrids and of introgressive hybridization events between two or more 
taxa (Hayashi et al., 2005), and of gene flow between populations (Alves et al., 2006; Ballard & 
Whitlock, 2004; Sánchez-Guillén, Wellenreuther, Cordero-Rivera, & Hansson, 2011). Moreover, 
phylogeographic studies of odonates have increased potential interest in the use of nuclear DNA 
sequence markers, as nuclear markers display a higher average diversity in invertebrates than in 
vertebrates (Bazin, Glémin, & Galtier, 2006). 
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Although the Coenagrionidae is the most diverse zygopteran family, surprisingly little is 
known about phylogeographic patterns of coenagrionids (but see Bernard et al., 2011; Jordan 
et al., 2005; Swaegers et al., 2014). Also, the group phylogeny remains a source of debate (Carle 
et al., 2008; Dijkstra et al., 2014; Dumont, Vierstraete, & Vanfleteren, 2010), highlighting the 
need for the development of new markers. Using these new nuclear markers in C. mercuriale, 
we uncovered high levels of polymorphism and considerable divergence among two geographic 
groups, a diversity that was not detected using 28S and H3. Polymorphism levels were higher 
in markers with greater length of intron in the fragment, underlining the perceived usefulness of 
sequencing introns for phylogeographic studies. On the other hand, as the primers are designed 
in exon regions (more conserved than introns) they can be potentially applied across a wide 
taxonomic range. 

The cross-amplification capacity of primers specifically designed for a species, in this case 
C. mercuriale, is expected to decrease with increasing genetic distance due to mutations in the 
primer sites, and is connected with the level of conservation of the exon fragment. This may 
explain why we could not obtain fragments of AgT and PGI-1219 for C. scitulum, while primers 
for PRMT amplified fragments from the genomes of zygopterans from three other families and 
even worked in species from a highly divergent group of gomphids (Table 3). Another possible 
explanation is that the targeted fragments of AgT and PGI-1219 in C. scitulum may be too long 
due to the presence of insertion(s) causing failure of amplification. Nevertheless, and given that 
more organisms have their genomes sequenced, the design of primers for these markers in other 
Odonata species would circumvent many of these potential limitations, providing a powerful and 
accessible tool for future studies. 

Moreover, Odonata phylogenetics and systematics studies would also benefit from an increase 
in the number of nuclear markers, as presently the markers available have recognized limited 
resolution power and many phylogenetic relations await clarification (Ballare & Ware, 2011; 
Blanke, Greve, Mokso, Beckman, & Misof, 2013; Dijkstra et al., 2014; Kim et al., 2014). 

In summary, we have developed five new nuclear polymorphic loci that can be used to improve 
our understanding of phylogeographic patterns of many coenagrionids in future studies, and 
one new locus that works on a wide taxonomic range of odonates. Odonata phylogeographic 
studies need to use several independent polymorphic nuclear DNA markers to provide sufficient 
resolution and statistical robustness to make novel insights into the evolutionary history and 
conservation genetics studies of dragonflies and damselflies. 
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